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Abstract 

Data collected at LEP at centre-of-mass energies = 189 — 209 GeV are used 
to study correlations of the spin of W bosons using e"^e~ W+W~ — iuqq events. 
Spin correlations are favoured by data, and found to agree with the Standard Model 
predictions. In addition, correlations between the W-boson decay planes are studied 
in e+e~ W"*"W~ — > iuqq and e"'"e~ — W"'"W~ — ^ qqqq events. Decay-plane 
correlations are measured to be consistent with the Standard Model predictions. 
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Introduction 



The study of the properties of the W boson constitutes one of the main physics goals of the 
LEP experiments. We have previously reported on the measurement of the helicity fractions 
of W bosons in e"'"e~ — ^ W^W~ — >■ iuqq events, with i denoting either an electron or a muon, 
using data collected at LEP at centre-of-mass energies, ^/s, up to 209 GeV [1]. Measurements of 
W-boson polarisation were also reported in the framework of a spin-density matrix analysis [2]. 
The same data sample is used here to study correlations of the spin of W bosons as well as 
correlations between the decay planes of the two W bosons. This last study also includes 
e"'"e~ — > W"'"W~ —>■ qqqq events. Some correlations of the spin of W bosons in the e"''e~ — > 
W^W~ process are expected in the Standard Model and are accessible with the size of the data 
sample collected at LEP. Weak decay plane correlations are also expected. 

Studies of triple-gauge-couplings of the W boson [3], which rely on the analysis of the full 
differential cross section of the e"'"e~ — >■ W"'"W~ process, have an implicit sensitivity to these 
correlations. Their results are presented in terms of couplings describing the Lorentz-invariant 
7WW and ZWW vertices. The aim of the analysis presented in this Article, as a natural 
extension of the work of Reference 1, is a direct and model independent measurement of the 
spin and decay-plane correlations. This enables a specific test of the corresponding Standard 
Model predictions, as suggested in References 4 and 5 for the spin and decay-plane correlations, 
respectively. Large deviations from the Standard Model predictions, possibly outside the reach 
of LEP, would suggest anomalous mechanisms for W-boson pair-production in e"''e~ collisions 
such as one-loop contributions from new particles or in general effects from any model with a 
symmetry breaking mechanism which affects interactions of longitudinally-polarised W bosons. 
To date, no such direct measurement was reported and this study introduces a new technique 
which could be exploited at future high-luminosity lepton colliders. 

W-boson spin correlations are measured in e"'"e~ — W"'"W~ iuqq events by tagging the 
helicity of the W boson which decays into hadrons and measuring the helicity of the W boson 
which decays into leptons. We consider two subsamples of W-boson pairs: the first is enriched 
in events where W~ bosons decaying into hadrons have a helicity Aw- = ±1 and the second is 
depleted of these events ^\ A difference in the helicity compositions of the W bosons decaying 
into leptons for the two subsamples would indicate the presence of W-boson spin correlations. 

The helicity fractions of the W~ boson decaying into leptons are obtained in a model inde- 
pendent approach from the shape of the distribution of the polar decay angle, 9}, between the 
charged lepton and the W~-boson flight direction in the W~-boson rest frame. The differential 
distribution of leptonic W~-boson decays is: 

^ /4 (l + cos^|)2 + /4 (l-cos^|)2 + /o^sin2^;, (1) 



iVdcos^l '^8 ' ■ K tj • 

where /_, /+ and /o represent the fractions of W~ bosons in the — 1, -|-1 and helicity states, 
respectively. Assuming CP invariance, these equal the fractions of the corresponding helicity 
states —1 and of the W"*" boson. For hadronic W-boson decays, considering only the 
absolute value of the cosine of the polar decay angle, | cos6'*|, the differential distribution is: 

^ '^^ /±^(l + |cos^*P) + /o^(l-|cos^*r), (2) 



A^d|cos^*| '^4' ' ^' ' •'^2' ' 



^^The charge conjugate state W+ is also included throughout this Article. CP conservation is assumed, as 
verified in W-boson polarisation studies [1,2]. 
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with f±=f^+f^. It is worthwhile to remark that the only hypothesis in the derivation of 
Equations 1 and 2 is the decay of a unit-spin W boson into fermions. 

The values of /+ and /o are a function of ^/s and the W~-boson production angle with 
respect to the incoming e~ , cos Qw- ■ The numerical values corresponding to the data sample 
investigated in this Article are therefore dependent on the integrated luminosity collected at 
each centre-of-mass energy at which the LEP collider was operated. 

The fractions of helicity combinations (Aw- , Aw+) depend strongly on cos Ow- • As an exam- 
ple, FigureC]shows the results of a leading-order analytical calculation [6]. This property is used 
to select intervals in which a particular helicity combination is enriched and thus spin correla- 
tions are more significant. Two intervals are considered: the forward bin, 0.3 < cos6w- < 0.9, 
where the fraction of the helicity combination (Aw-,Aw+) = ("l^+l) is increased to about 
68% of all W-boson pairs, compared to an average value of 49% over the whole cos 6w- range; 
and the backward bin, —0.9 < cos6w- < —0.3, where the fraction of the helicity combination 
(Aw-, Aw+) = (0, 0) is increased to about 26%, compared to an average value of 8%. 

To tag the helicity of the W boson which decays into hadrons, cuts on | cos^*| are used. 
According to Equation 2, for small values of | cos 6**1, the sample is depleted of helicity ibl 
states, while for large values of | cos^*| the sample is enriched in helicity ±1 states, as shown 
in Figure El The interval |cos^*| < 0.33 is chosen for the Aw- = ±1 depleted sample and 
the interval |cos6'*| > 0.66 for the Aw- = ±1 enriched sample. The helicity fractions of the 
W~ bosons decaying into leptons are obtained from a fit of Equation 1 to the distribution of 
cos^; [7]. 

Decay-plane correlations are studied in both e"'"e~ — > W^W" iuqq and e''"e^ W''"W^ ^ 
qqqq events using the absolute value of the angle, |A0|, between the planes defined by the de- 
cay products of the two W bosons in the rest frame of the W-boson pair. The strength of the 
correlation is measured by the parameter D of the differential distribution [8]: 

The experimental distribution of |A0| is fitted to obtain the parameter D [7]. The expected 
value of D in the Standard Model, for the ^/s range under investigation, is estimated with 
a leading-order analytical calculation [8] to vary from 0.021 at ^/s = 189 GeV to 0.015 at 
■\/s = 209 GeV, with a luminosity-weighted average of 0.018. 

For both analyses, prior to the fits, the distributions of the measured angles are corrected 
by means of an efficiency correction function which accounts for selection efficiencies, migration 
effects and the presence of initial state radiation. This correction function is obtained from 
large Monte Carlo samples as the ratio of the angular distributions for simulated and generated 
events. 

Data and Monte Carlo Samples 

The analysis is based on 629.3 pb~^ of data collected with the L3 detector [9] at ^/s = 189 — 
209 GeV, as detailed in Table ^ corresponding to a luminosity-weighted average y/s value of 
197.9 GeV. 

Signal Monte Carlo events are generated using the KORALW [10] program for the e~^e~ — >■ 
W^W~ — s> euqq, /ii/qq and qqqq processes. The Standard Model predictions for the W-boson 
spin correlations and the decay-plane correlations are obtained from the distributions generated 
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at different values of ^/s, combined according to the collected luminosity. Background processes 
are generated using KORALW for all other final states of W-boson pair-production, KK2f [11] 
for the e"'"e~ qq(7) process, and PYTHIA [12] for the e"''e~ ZZ process. For studies of 
systematic effects, signal events are also generated using the EEWW [13], YFSWW [14] and 
EXCALIBUR [15] programs. 

To test the consistency of the predictions for W-boson spin and decay-plane correlations, 
large samples of signal events are generated using the EEWW and YFSWW Monte Carlo 
programs. These differ from KORALW in the level and implementation of 0{a) corrections. 
The predictions given by the three programs for the strength of the correlations are in agreement 
with each other, within their own statistical uncertainties. 

The predicted Standard Model value of the parameter D, which describes the decay-plane 
correlations is D = 0.010±0.002, as obtained with the KORALW Monte Carlo. The uncertainty 
is statistical. This differs from the value D = 0.018 obtained by an analytical calculation [8] 
due to the inclusion of radiative effects in the Monte Carlo. 

The L3 detector response is simulated with the GEANT [16] and GHEISHA [17] packages. 
Detector inefficiencies, as monitored during the data taking periods, are included. 

After detector simulation, the average resolution on cosBw- is found to be 0.06 while the 
resolution on cos 6'* is found to be 0.10 for electrons, 0.11 for muons and 0.14 for hadrons. The 
resolution on |A0| is 7.5° for semi-leptonic events and 10.0° for hadronic events. 

Event Selection 

The selection of e"'"e~ — W~'~W~ — iuqq events is the same as for the study of the W- 
boson polarisation [1], however, the low-statistics data sample collected at = 183 GeV 
is not considered here. The numbers of selected events are listed in Table ^ for different 
values of y/s. In total, 1861 events are selected with an average efficiency of 67.6% and an 
average purity of 96.6%, which are only slightly dependent on ^/s. The contamination from 
the e+e" — > W"'"W~ — > Tuqq and e"'"e~ qq(7) processes is 2.3% and 1.1%, respectively. 

The selection of e"'"e~ — > W'''W~ qqqq events is performed as follows. High multiplicity 
events are selected by requiring more than 20 charged tracks and more than 25 calorimetric 
clusters. The visible energy of the event must satisfy E^is > 0.75^/s. To reject the e"'"e~ — > qq(7) 
background, the event thrust must be less than 0.88 and the polar angle of the thrust axis, 
6t, has to satisfy |cos^t| < 0.95. The missing momentum of the event has to be less than 
60 GeV. Events containing electrons, muons or photons with energy greater than 20, 20 or 40 
GeV, respectively, are rejected. Jets are reconstructed using the Durham algorithm [18], with a 
jet-resolution parameter for which the event changes from a four-jet into a three-jet topology, 
1/34, greater than 0.0015. Two pairs of jets are formed, corresponding to two W bosons. Of 
the three combinations, the optimal pairing of jets is chosen as the one with the smallest mass 
difference, disregarding the pairing corresponding to the smallest mass sum. This algorithm 
yields the correct assignment of jets to W bosons for about 70% of the selected Monte Carlo 
events. Finally, the reconstructed W bosons must have a mass between 40 and 120 GeV. 
Figure |H1 shows the distributions of some selection variables for data and Monte Carlo. 

The numbers of events selected by these criteria at different values of ^/s are listed in 
Table lU In total, 4919 events are selected with an average efficiency of 76.3% and an average 
purity of 75.7%, which do not strongly depend on ^/s. The background contamination is 18.9% 
from the e"'"e~ — >■ qq(7) process, 4.9% from the e"'"e~ ZZ process and 0.5% from W-boson 
pairs which decay into other final states. 
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Analysis of W-Boson Spin Correlations 



The values of cos^| and |cos^*| are determined for each selected event. The latter is ap- 
proximated by the absolute value of the cosine of the angle of the thrust axis of the W-boson 
decaying into hadrons, calculated with respect to the W-boson flight direction in the W-boson 
rest frame. 

The W-pair events are then classified according to the values of cosGw- and |cos^*| to 
build four test samples: 

O.3<cos0w-< 0.9, I cos^*| < 0.33 ; 0.3<cosew-< 0.9, | cos > 0.66 ; 

-0.9 < cosGw- < -0.3, I cos^*| < 0.33 ; -0.9 < cosGw- < -0.3, | cos^*| > 0.66 . 

In each of the samples, the fractions of W-boson helicity states for W bosons decaying into 
leptons are obtained from the event distributions, dA^/dcos6'^. For each energy point the back- 
ground, as obtained from Monte Carlo simulations, is subtracted from the data. The corrected 
decay angle distributions at different values of ^/s are combined into single distributions, shown 
in Figure I3J which are then fitted with the function in Equation 1, using /_ and /o as the fit 
parameters. The fraction /+ is obtained by constraining the sum of all three parameters to 
unity. 

Finally, the fitted fractions are corrected for the bias which originates from migration effects 
due to detector resolution [1, 19]. Bias correction functions are determined in each of the 
investigated cosGw- and | cos6'*| ranges. They vary from 1% to 15%. 

Several sources of systematic uncertainty are considered, as summarised in Table 121 for the 
measurement of the helicity fractions /_ and /o for all bins of cos 0w- and | cos 6** | . All selection 
cuts are varied over a range of one standard deviation of the corresponding reconstruction 
accuracy. The corresponding variation of the helicity fractions, corrected for its statistical 
component, is taken as systematics. All fits are repeated with one bin more or one bin less in 
the angular distributions. The average difference is retained as systematics. Uncertainties on 
the bias and the efficiency corrections are determined by varying the bias correction function 
and the efficiencies in each bin by one standard deviation, as derived from the statistics of the 
corresponding Monte Carlo samples. Additional contamination from the non double-resonant 
four-fermion final states is evaluated using the EXCALIBUR Monte Carlo. Background levels 
are varied according to Monte Carlo statistics for the background processes. 

The results of the fits are summarised in Table El The correlation coefficients of the param- 
eters /_ and /o derived from the fit are shown in Table IH W-boson spin correlations would 
manifest as sizable differences between the values of /_ and /o measured for the samples de- 
pleted of and enriched in ±1 helicity, respectively. These differences are also listed in Table El 
There is a general agreement between the observations and the expectation. The largest differ- 
ence is observed for /_ in the forward bin, with an observation of —0.32 ± 0.12 to be compared 
with a prediction of —0.11 ± 0.01. Two consistency tests are performed. In the first, only data 
is considered and a confidence level is calculated for the absence of W-boson spin correlations. 
This confidence level is 1.2%, which allows to conclude that W-boson spin correlations are 
observed. A second test compares the data with the Standard Model KORALW Monte Carlo. 
The confidence level for their compatibility is 34.7%. 
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Analysis of Decay-Plane Correlations 



Events from the e+e~ — W"'"W~ — > iuqq and e^e~ — W"'"W~ — qqqq processes are used to 
study correlations between the W-boson decay planes. 

For e"'"e~ — >■ W"'"W~ — > iuqq events, the neutrino momentum is derived from the total 
missing momentum of the event. The decay plane of the W boson decaying into leptons is 
determined from the lepton and the neutrino directions. The decay plane of the W boson 
decaying into hadrons is determined from its thrust axis in the W-boson rest frame and the W- 
boson flight direction. For e+e~ — > W+W~ — > qqqq events the reconstructed jets are boosted 
into the W-pair rest frame and the decay planes are determined by the two jets assigned to 
each reconstructed W boson. The angle |A0| between the decay planes of the two W bosons 
is then calculated. For each value of ^/s, the |A0| distribution is corrected for efficiency and 
background, taking also into account wrongly-paired four-jet events. The efficiency corrections 
are found to be basically independent from the values of D. The corrected distributions are 
combined into a single distribution, shown in Figure El A binned fit for D, using Equation 3, 
is performed on the normalised distribution. 

The same systematic studies are performed as for the W-boson spin correlations. Addition- 
ally, for e"''e~ — > W+W~ qqqq events, several pairing algorithms are used as a cross-check. 
The largest difference in the fit result between the pairing methods is taken as systematic un- 
certainty. To reproduce the measured four jet event rate of the e"'"e~ — ^ qq(7) background [20], 
the corresponding Monte Carlo is scaled by -|-10%. Half of the effect is retained as systemat- 
ics. Table El summarises the systematic uncertainties in the measurement of the decay-plane 
correlation parameter D. 

The resulting value of D for e'^e^ W^W~ — > iuqq events is found to be 0.051 ± 0.033 ± 
0.019, where the first uncertainty is statistical and the second systematic. It is in agreement 
with the Standard Model prediction from the KORALW Monte Carlo of = 0.006 ± 0.004, 
where the error refiects the Monte Carlo statistics. For e"'"e~ — > W^W~ — > qqqq events, 
D is found to be —0.016 ± 0.028 ± 0.016, in agreement with the KORALW prediction of 
D = 0.013 ± 0.003. Combining the two decay channels, a value D = 0.012 ± 0.021 ± 0.012 
is found in data, in agreement with the combined value from the Standard Model KORALW 
Monte Carlo oi D = 0.010 ± 0.002. 

Summary 

In conclusion, this study completes our investigation of W-boson polarisation. A new technique, 
which is promising for future high-luminosity electron-positron colliders, is deployed and allows 
to establish the existence of W-boson spin correlations. Their magnitude is found to agree with 
the Standard Model predictions. In addition, W-boson decay-plane correlations are studied for 
the first time, and no large deviations with respect to the Standard Model predictions, which 
could hint to New Physics, are observed. 
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Table 1: Average centre-of-mass energies with corresponding integrated luminosities 
and numbers of selected events. 



0.3 < cosOw- < 0.9 



±1 depleted 







±1 enriched 



/o 



-0.9 < cosGw- < -0.3 



±1 depleted 







±1 enriched 







Selection 
Fit binning 
Bias correction 
Efficiency correction 
Four fermions 
Background levels 



0.034 
0.027 
0.018 
0.003 
0.001 
0.006 



0.045 
0.021 
0.026 
0.001 
0.006 
0.008 



0.010 
0.030 
0.010 
0.005 
0.006 
0.004 



< 0.001 
0.057 
0.013 
0.004 
0.003 
0.008 



0.053 
0.078 
0.031 
0.004 
0.012 
0.022 



0.048 
0.097 
0.055 
0.006 
0.024 
0.024 



0.050 
0.020 
0.016 
0.002 
0.003 
0.014 



0.087 
0.051 
0.034 
0.005 
0.026 
0.019 



Total 



0.048 



0.057 



0.034 



0.059 



0.102 



0.126 



0.058 



0.111 



Table 2: Systematic uncertainties on the measurement of /_ and /o for leptonic 
W-boson decays in bins of cos Ow- cind for samples depleted of and enriched in the 
±1 helicity. The statistical component of these uncertainties is removed. 
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0.3 < cos ©w^ 


< 0.9: fAw~, Aw+) = (- 


1,±1) enriched 




f 

J — 


J + 


fn 
J u 


Data ±1 depleted 


0.521 ±0.086 ±0.048 


0.121 ±0.058 ±0.035 


0.358 ±0.116 ±0.057 


Data ibl enriched 


0.839 ±0.057 ±0.034 


0.087 ±0.034 ±0.036 


0.074 ±0.072 ±0.059 


Difference 


-0.318 ± 0.103 ± 0.059 


0.034 ± 0.067 ± 0.050 


0.284 ±0.137 ±0.082 


MC ±1 depleted 


0.670 ±0.008 


0.148 ±0.006 


0.182 ±0.012 


MC ±1 enriched 


0.781 ±0.007 


0.091 ±0.004 


0.128 ±0.009 


Difference 


-0.111 ±0.011 


0.057 ±0.007 


0.054 ±0.015 




—0.9 < cosOw 


— *\ — (j.o! ( \\\r— . \\AT+ ) = 


- (0, 0) enriched 




f 


J + 


fn 

J u 


Data ±1 depleted 


0.387 ±0.136 ±0.102 


0.398 ±0.121 ±0.061 


0.215 ±0.316 ±0.126 


Data ihl enriched 


0.152 ± 0.082 ± 0.058 


0.530 ±0.101 ±0.087 


0.318 ± 0.215 ± 0.111 


Difference 


0.235 ± 0.159 ± 0.117 


-0.132 ± 0.158 ± 0.107 


-0.103 ± 0.382 ± 0.168 


MC ±1 depleted 


0.186 ±0.014 


0.351 ±0.016 


0.463 ± 0.023 


MC ±1 enriched 


0.158 ±0.014 


0.532 ±0.017 


0.310 ±0.026 


Difference 


0.028 ±0.020 


-0.181 ±0.023 


0.153 ±0.035 



Table 3: The W-boson helicity fractions measured for different intervals of cosGw-- 
The results are shown for different subsamples depleted of or enriched in helicity ±1. 
The corresponding helicity fractions in the Standard Model, as implemented in the 
KORALW Monte Carlo program, are also given with their statistical uncertainties. 





/_ — /o correlation coefficient 


Data 


Monte Carlo 


±1 depleted 


±1 enriched 


±1 depleted 


±1 enriched 


0.3 < cosGw- < 0.9 


-92% 


-93% 


-88% 


-91% 


-0.9 < cosOw- < -0.3 


-76% 


-77% 


-80% 


-78% 



Table 4: Correlation coefficients between the fit parameters /_ and /q. 





e+e- ^ W+W- - 




e'^e - 


W+W- - 


qqqq 


Selection 


0.013 






0.007 




Fit binning 


0.012 






0.007 




Efficiency correction 


0.001 






0.002 




Four fermions 


0.005 










Background levels 


0.002 






0.007 




Jet pairing 








0.011 




Total 


0.019 


0.016 



Table 5: Systematic uncertainties on the measurement of the decay-plane correlation 
parameter D. The statistical component of these uncertainties is removed. 
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Figure 1: Differential cross section for pair production of polarised W bosons at ^/s = 190 
GeV averaged over initial electron polarisations. The W~ and W"*" helicities in the e"'"e~ centre- 
of-mass frame are given in parentheses. The intervals used for the analysis of W-boson spin 
correlations are indicated by the arrows. 
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Figure 2: Standard Model predictions for the relative contribution of the helicity states ±1 to 
the e^e~ — > W"'"W~ — iuqq differential cross section as a function of | cos^*|. The intervals 
used for the analysis of W-boson spin correlations are indicated. 
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Figure 3: Distributions of variables used for the selection of e"'"e~ W^W~ 
qqqq events: (a) normalised visible energy, (b) event thrust, (c) energy of the most 
energetic photon, E^, (d) jet-resolution parameter, — log 2/34. In each plot, all other 
selection criteria are applied. The arrows indicate the positions of the cuts. 
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Figure 4: Corrected cos 9} distributions for W £z/ decays for data and the KORALW Monte 
Carlo in the intervals (a) 0.3 < cosG-^- < 0.9, and (b) —0.9 < cos6-\y- < —0.3. The 
distributions are shown after classifying the events in two samples, according to the helicity of 
the W boson decaying into hadrons. The first sample is enriched in W bosons in a transverse- 
helicity state, the second sample is depleted of W bosons in a transverse-helicity state. For 
clarity, the data points are slightly shifted. 
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Figure 5: Corrected |A0| distributions for (a) e+e~ W"^W~ — » iuqq and (b) e+e~ — >• 
W"''W~ — > qqqq events for data and the KORALW Monte Carlo. The fit results are also 
shown. 
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